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Abstract — In this paper, we investigate two new candidate 
transmission schemes, Non-Orthogonal Frequency Reuse (NOFR) 
and Beam-Hoping (BH). They operate in different domains 
(frequency and time/space, respectively), and we want to know 
which domain shows overall best performance. We propose 
a novel formulation of the Signal-to-Interference plus Noise 
Ratio (SINR) which allows us to prove the frequency/time 
duality of these schemes. Further, we propose two novel capacity 
optimization approaches assuming per-beam SINR constraints in 
order to use the satellite resources (e.g. power and bandwidth) 
more efficiently. Moreover, we develop a general methodology to 
include technological constraints due to realistic implementations, 
and obtain the main factors that prevent the two technologies 
dual of each other in practice, and formulate the technological 
gap between them. The Shannon capacity (upper bound) and 
current state-of-the-art coding and modulations are analyzed in 
order to quantify the gap and to evaluate the performance of 
the two candidate schemes. Simulation results show significant 
improvements in terms of power gain, spectral efficiency and 
traffic matching ratio when comparing with conventional systems, 
which are designed based on uniform bandwidth and power 
allocation. The results also show that BH system turns out to 
show a less complex design and performs better than NOFR 
system specially for non-real time services. 

Index Terms — Multibeam Satellite, Duality, Beam-Hopping, 
Frequency-reuse, and Time-reuse. 

I. Introduction 

THE efficient management of multibeam antenna re- 
sources, e.g. power, bandwidth and time-slot, is crucial 
for economic competitiveness. Specifically, in modern satellite 
networks, each satellite uses multiple beams, each of which 
illuminates a cell on the ground to serve a coverage area. 
Multibeam antenna technology is used because it can increase 
the total system capacity significantly [l]-[3]. However, each 
beam will compete with others for resources to achieve 
satisfactory communication. This is due to the fact that the 
traffic demand throughout the coverage is potentially highly 
asymmetrical among the beams. Therefore, the satellite re- 
quires a certain degree of flexibility in allocating the power, 
bandwidth and time-slot resources to achieve a good match 
between offered and requested traffic. 

There are some precedents of resource allocation optimiza- 
tion techniques for satellite systems. In [4], [5], the authors 
investigate the dynamic bandwidth allocation techniques, and 
in [6] the authors propose a quasi-optimal solution to manage 
the frequency slots allocation to service providers, however, 
the results are only for the satellite uplink. A power allocation 
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policy is proposed in [7], which suggests to stabilize the 
system based on the amount of unfinished work in the queue 
and the channel state, and a routing decision is made for 
the maximum total throughput. However, the authors do not 
taking into account the co-channel interference. In [8], [9], a 
tradeoff strategy is proposed between different objectives and 
system optimization. The power allocation is optimized based 
on the traffic distribution and channel conditions. However, 
the co-channel interference is not taken into account. In [10], 
a joint power and carrier allocation problem is discussed, 
however, it focus on the return (RTN) uplink. In [11], [12], 
the authors focused on the capacity optimization in multibeam 
satellite system, and the duality in frequency and time domain 
is studied. The optimization problem of power and carrier 
allocation has also been addressed in terrestrial networks. E.g., 
the authors in [13] propose an axiomatic-based interference 
model for Signal-to-Interference plus Noise (SINR) balancing 
problem, but the conclusions cannot be directly extrapolated 
to a satellite scenario. Although the resource allocation op- 
timization has been study extensively, the objective of this 
paper is different from the aforementioned literatures in var- 
ious aspects. E.g., most of the existing resource allocation 
optimization approaches focus on terrestrial networks or on the 
satellite RTN link, while we focus on the satellite FWD link. 
In addition, we address the resource allocation according to the 
realistic asymmetric traffic distribution by managing the co- 
channel interference due to frequency reuse. Existing work has 
focused on the analysis of the resources allocation in frequency 
domain. Our aim is to characterize the best resource allocation 
scheme in multi-domain, and to show in which domain the 
overall performance is best. 

In this paper, we investigate two new transmission schemes, 
Non-Orthogonal Frequency Reuse (NOFR) technique and 
Beam-Hopping (BH), which have been chosen as candidates to 
replace current regular frequency reuse transmission scheme. 
The first one is designed based on the frequency division over 
a flexible payload design which allows managing interference 
as an alternative to a complete orthogonal frequency reuse. 
The second one is based on the time/space division. Both 
techniques can potentially cope with the asymmetric traffic 
distribution as opposed to current satellite resources allocation 
scheme, which is designed to allocate fixed power and band- 
width to each ground cell. This leads to a waste of resources 
in low traffic requirement beams. On the contrary, it does not 
satisfy traffic demand in the hot ground cells, where the traffic 
requirement is high. 

In this paper, we study the duality between NOFR and 
BH, i.e. frequency and time/space duality. The concept of 
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duality gives rise to many interesting properties to simplify 
system models. Generally speaking, a duality translates con- 
cepts, theorems or mathematical structures into other concepts, 
theorems or structures, in a one-to-one fashion, often (but 
not always) by means of an involution operation: if the dual 
of A is B, then the dual of B is A. E.g., the duality of 
space and time is studied in [14], [15], Gaussian multiple- 
access/broadcast channel duality is discussed in [16], uplink 
and downlink duality is presented in [17], [18]. We develop a 
general methodology to study the duality of the two schemes 
that also considers the technological constraints due to realistic 
implementations, and obtain the main factors that prevent the 
two schemes be in practice dual of each other. 

The novel contributions of this paper can be summarized as 
follows: 

* The frequency and time duality is formulated for the 
multibeam satellite system, and the duality conditions are 
derived for a practical system. 

• We prove that new transmission schemes, NOFR and BH, 
can match much better than the conventional design in the 
realistic asymmetric traffic model, and also prove that the 
BH system performs only slightly better than NOFR. 

The rest of this paper is organized as follows: In Section II, 
the problem statement is presented. In Section III, we model 
the multibeam downlink system to obtain a mathematical 
expression of SINR. The duality of NOFR and BH is discussed 
in Section IV. In Section V, we formulate and solve the 
satellite capacity optimization problems. In Section VI, the 
technological gap is obtained with a realistic system payload 
model. The simulation results are presented in Section VII. In 
Section VIII, we draw the conclusions of the paper. 

II. Problem Statement 

In multibeam satellite systems, the beamforming antenna 
generates K beams over the coverage area. For both NOFR 
and BH systems, we firstly introduce some payload parameters 
(as shown in Table I). 

« Granularity: B c is the carrier granularity defined as B c = 
B tot /N c in NOFR system. It means that the allocated 
bandwidth per ground cell should be an integral multiple 
of B c . We use T s , with the same meaning but in BH 
system, i.e. the minimum unit of time duration that can 
be allocated per cell. 
> Resource allocation matrix: Wij and tij are the elements 
of the resource allocation matrix for the NOFR and BH 
systems, respectively. The matrix indicates which carrier 
or time-slot j is allocated to the ground cell i. Note that 
BH can direct the satellite beams to specific ground cells, 
i.e. it is a space allocation too. 
In the case of a NOFR scheme, each ground cell can be 
allocated a variable number of carriers (e.g. Ni, as shown 
in Fig.l) depending on the traffic requirement. Carriers can 
be re-used throughout the coverage, but we do not impose 
any restrictions on the frequency reuse, it will be given by 
the resource optimization (i.e. interference minimization for 
a given traffic demand pattern) and therefore will be non- 
orthogonal. In the case of a BH system, the total bandwidth is 
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Fig. 1. Bandwidth segmentation. 

simultaneously used by a set of the ground cells during a time- 
slot (T s ). We assume that the resource allocation takes place 
during a given time window divided into N t time-slot. Each 
ground cell can be allocated a variable number of time-slot. 

Note that both techniques allow a number of ground cells 
to use the same frequency band or time-slot, resulting in co- 
channel interference. The problem tackled in this paper is to 
optimize the capacity by taking into account the co-channel 
interference. Further, we prove the duality of these techniques 
by developing a formulation that also allows including techno- 
logical constraints. Moreover, we compare the performance of 
the proposed new transmission schemes with the current one 
for the realistic asymmetric traffic model. To do so, we propose 
a novel iterative algorithm, which do not only optimize the 
power and bandwidth allocation (for NOFR systems), but 
also optimize the structure of the spectral mask matrix W. 
This matrix indicates which carriers are allocated per-beam in 
order to minimize the co-channel interference. Although the 
power and carrier optimization problem has been addressed in 
terrestrial networks, it is new in satellite communications. 

III. Multibeam System Model 

In this section, we first formulate the multibeam system 
model in frequency domain (i.e. for the NOFR scheme). After 
that we state the conditions for duality and prove that NOFR 
and BH are dual of each other and hence, the formulation is 
also valid in time domain(i.e. for the BH scheme). This dual 
formulation allows us to derive a unique SINR expression, 
which will be used in the following section for capacity opti- 
mization. Following, we introduce the different sub-models. 

A. Channel Model 

We do an analysis in time and hence the channel attenuation 
corresponds to the free space losses and atmospheric losses (in 
case of frequencies above Ka band). We assume an instanta- 
neous analysis with fixed coefficient. The channel attenuation 
amplitude matrix A G C A xA is defined as 

A = diag{ai,a 2 , ■ • ■ ,a K } , (1) 

where ai denotes the channel attenuation factor over the 
destination user beam i. 

B. Antenna Model 

An Array Feed Reflector (AFR) based Antenna system is 
assumed, it can generate a regular beam grid array consisting 
of a very high number of highly overlapping, narrow beam 
width, composite user beams. Each beam is synthesized by 
adding array elements whose phases and amplitudes are ad- 
justable, and hence we can provide flexible power allocation 
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by controlling the On-Board Processor (OBP). Therefore, we 
suppose that the antenna gain matrix G E 



KxK is given as 



flu 
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(2) 



where \g^\ 



:1 is the antenna gain of the on-board 



antenna feeds for the jth beam towards the ith user beam. 
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Frequency-time duality 
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C. Received Signal Model 

In the frequency domain, the transmitted symbols over 
N c carriers to beam i (i = 1,2, ••• ,K) are defined as 
Xj = [xn , Xi2 , • • ■ ,XiArJ T . Let the spectral mask matrix 
W E R N ° xK be defined as W = [wi,w 2 ,--- ,w K ], and 
the ith column vector w» E R 7VcXl be defined as w, = 
[wn,Wi2, ■ ■ ■ , WiN a ] , which is the spectral mask vector for 
beam i and indicates which TDM carriers and how much 
power is allocated to beam i. 

Let H = AG be the overall channel matrix, and W t = 
diag {wi}. Then the received signal by all the iV c carriers for 
beam i, y, E C JV ° xl , can be expressed as desired signal and 
interference as 

K 

y t = h vl x t + h ik*k + n 4 , (3) 

fc = l(k^i) 

where Xj is the spectral masked symbol vector for beam i, 
defined as xj = WjX,. The first term corresponds to the 
desired signals coming from the ith on-board antenna. The 
second term is the sum of interference signals from the other 
on-board antennas, E C^* 1 is a column vector of zero- 
mean complex circular Gaussian noise with variance a 2 at 
beam i. 

D. Signal-to-Interference plus Noise Ratio 

In the frequency domain, the bandwidth is segmented as 
shown in Fig.l. We assume that the whole bandwidth is 
segmented in N c carriers. The spectral mask matrix can be 
reformulated as W = [w^wj,--- ,w^] T , where Wj = 
[wij, u>2j, ■ ■ • , wxj], indicates which beams are allocated 
carrier j. Let the ith row of H be defined as = 
[hii,h i2 ,--- ,h iK ] and hi = h.i\ (hu = 0) is the channel of 
interference contribution. We assume that the amplitude of 
the transmitted symbols is normalized (i.e. \xij\ 2 = 1, Vi = 
A:.'./ 1.---..V.). 

Then, the transmitted signal power of all the carriers for 
beam i can be given by the diagonal elements of the matrix 
S{ E M. N " xNc as (note that the superscript / and t indicate 
the expression in frequency and time domain, respectively) 

Sf = IfciifWiW?. (4) 

And the co-channel interference power of all the carriers 
for beam i can also be given by the diagonal elements of the 
matrix Uf £ R N ^ N - as 

Uf = diag { [hiwf w,hf ] . =i 2i ^ } . (5) 



Thus, the interference power plus the noise matrix, , will 
be given as 

Vf = Uf +a 2 I Nc . (6) 

Consequently, the SINR for beam i, defined as r{ E 

M. NcXN,: , can be expressed as 

rf = sf(vf)- 1 . (?) 

Obviously, T{ is a diagonal matrix, because both S{ and 
V{ are diagonal matrixes. Thus, the SINR for the jth carrier 
used by beam i will be the jth diagonal element of the matrix 
T{. This means that for each carrier j of beam i, the SINR 
can be formulated as 



^ \h lk w k3 \ 2 + a 2 
fc=i(fc^») 

For the RTN uplink scenario, the authors in [10] formulate 
the SINR in a similar way for a specific terminal. However, 
for the FWD downlink, we formulate the SINR per beam, and 
all the carriers' SINRs are integrated in an equation (7), the 
SINR for a specific carrier (8) is also derived from (7). 

IV. Frequency/Time Duality 

In the previous section, expression (8) gives the SINR in 
terms of the spectral mask vector (i.e., in frequency domain). 
In this section, we will propose the frequency/time duality of 
(8). For doing so, we first state the dual expression in time 
domain. After that, we can find the conditions for the duality. 

A. Dual System Model 

In the time domain, the time window is segmented into 
Nt time-slot. The time-slot mask matrix can be formulated as 
T = [tf ,t|\ • • • ,t^J T , where ij = [tij,t 2 j, ■ ■ ■ ,t K j], indi- 
cates which beams are allocated time-slot j. Then, according 
to the duality elements in Table I, the transmitted signal power 
matrix S*, the co-channel interference power matrix U*, the 
interference power plus the noise matrix V* and the SINR 
matrix T\ in time domain can be formulated as follows 

Sj = |M 2 T,Tf, (9) 
U* = diag{[Mft,hf] . =i2i (10) 
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Vj = U* + cr 2 I 



rj = sj(vj) 



(ID 



(12) 



The SINR of beam z and time-slot j will be the jth diagonal 
element of the matrix T\. Hence, the SINR can be formulated 

as 

| hatij | 



(13) 



E 



\hiktkj | 2 + o" 2 



fe=i(fc^») 

From a theoretical point of view [14], (8) and (13) are dual 
of each other. However, for a practical system, we derive the 
duality conditions in the next section. 

B. Duality Conditions 

From (8) and (13) we can extract the duality conditions. 
In order to do so, we first express the beam-level sum-rate 
throughput as follows 



and the dual is 



N t 

EBtot t 

3 = 1 



(14) 



(15) 



where jjy = f(jij) is the spectral efficiency (ry^ can be 77/- 
or ijlj in frequency domain or time domain, respectively), and 
/(7»j) equals to log 2 (l+7ij) for Shannon limit with Gaussian 
coding, or can be a quasi-linear function in DVB-S2 [19] with 
respect to SINR . 

Hence, in order to obtain the duality conditions, we assume 
that R{ = R\ and the throughput rate per-carrier in frequency 
domain (or per time-slot in time domain ) is equivalent for 
each illuminated beam. The following conditions should be 
fulfilled for systems to be dual in practice: 

• Granularity in frequency and time domains should be the 
same: N c — N t . 

• The entries of resource allocation matrix should be the 
same in frequency and time domains: uiij = tij. 

• The spectral efficiency function /(■) should be the same 
for NOFR and BH systems in frequency and time do- 
mains, respectively: f(r]{j) = f(r)ij)- 

V. Capacity Optimization 

In this section, we propose two capacity optimization prob- 
lems under the constraints of the traffic requested per-beam 
and overall available power. The first problem, PI: capacity 
optimizing with co-channel interference, maximizes the total 
capacity allocated with respect to the traffic requested. Since 
the problem is not only non-convex but also to the need 
of preserving the geometry of the matrix W. Therefore, we 
propose an iterative algorithm solution. The second one, P2: 
capacity optimizing without co-channel interference, is a sim- 
plified problem by assuming that the co-channel interference is 
negligible, hence, the problem can be solved straightforwardly 
by the lagrangian approach. 



TABLE II 

Algorithm solution for NOFR system 

1: Initialize: R k <= 0, VTc. n it <= 0. W <= 
2: i <= 0. 

Generating beam set A s : 

As = {ti,»a,-- - ,*jv|0 < 5^ < < l}. 

where i n € {1,2, ••• ,K},n'= 1,2," ■ , N. 
3: nit <= nn + 1 
Repeat: i i + 1. k <= A s {i) 
4: Solve the Rayleigh quotient problem: 



5: w kj 



are max — 1= — K-^- 

° e H v'e, 

Nl/2 



6: Update S f k , V f k <= Jj{ + a 2 I 

7: go to step 3, 
until k > in- 

8: Update j f kj ,Vk,j. R k <= £f=i fjij&.Vfc. 

9: go to step 2, until As is empty or ^^j, wfwj < P t 



A. PI: capacity optimizing with co-channel interference 

Obviously, 7^ in formula (8) not only depends on the 
spectral mask vector of beam i (Wj), but also depends on 
that of the co-channel beams. And hence, the spectral mask 
vector for each beam must be optimized jointly with the others. 
The specific design of one beam's spectral mask vector may 
affect the crosstalk experienced by other beams. Hence it's 
a complicated task to design the spectral mask matrix W 
jointly. In order to best match the offered and the requested 
traffic per-beam, we develop a methodology to solve the 
spectral mask matrix W in this section and to jointly optimize 
power and carrier allocation. Note that we only discuss the 
capacity optimization for NOFR system because BH is dual 
with NOFR, thus the formulation is also applicable for BH 
system by changing the duality parameters in Table I. 

Existing results in the references [20]-[23] is exclusively 
over the power allocation. However, we assume an additional 
degree of freedom: carrier allocation (bandwidth granularity). 
We propose to use binary power allocation (BPA) (|wij| 2 = 
{0,-P S at},* = 1, 2, • • • , K;j = 1,2, •■■ ,N C ) and quantized 
bandwidth allocation in order to decrease the complexity, 
where P sat is the TWTA saturation power per carrier. 

1) Optimization Problem Formulation: The capacity opti- 
mization problem can be formulated as 



max 
w 



K 

E 



Ri 



subject to Ri < Ri, 

K 

^wf Wj < P tot ; and 

»=i 



(16) 



w 



{0,Psat},Vi,j. 



where Ri is the traffic requested by beam i, Pj(W) is defined 
in Table I. P tot is total available satellite power, P sat is the 
saturation power per carrier, which is limited by the satellite 
amplifier. 

2) Iterative Algorithm Solution: The general analytical so- 
lution of (16) is a complex problem due not only to the clear 
non-convexity but also to the need of preserving the geometry 
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of the optimization model (i.e. the structure of matrix W). 
Therefore, we propose an iterative algorithm solution, which 
is summarized in Table II. The beam set A s is constituted by 
all the beams, in which the traffic requirement is not achieved 
(i.e. £fc < 1). Quantities associated with the nth iteration are 

rift 

denoted by n it . Each iteration is based on a two-step process. 

Firstly, we optimize subspace-by-subspace and obtain an 
analytical solution to the sub-problem of allocating the carrier 
on a per-beam basis (as shown in step 4 of Table II). The 
optimal carrier allocation per-beam can be formulated as a 
generalized Rayleigh quotient, e.g. for beam i, the problem 
can be formulated as: 

arg max — f — , 

3 efV/ e, 



K 



subject to wf^Wi < P to 



(17) 



where e.,- 6 



i Nc x 1 is standard basis vector, which denotes the 
vector with a "1" in the jth coordinate and O's elsewhere. 

The solution of the generalized Rayleigh quotient problem 
shown in (17) is given as 



«max(S{(V/)- 1 )= Wmax (r{), 



(18) 



where v m - dx (T{) (as expressed in 7) indicates the eigenvector 
related to the maximum eigenvalue of matrix Tj . 

Secondly, we obtain the power allocated to the selected 
carriers from the power constraint (as shown in step 5 of Table 
II). Wij for the jth carrier of beam k can be obtained with the 
solution of e. as 



H it, \l/2 



(19) 



After each iteration, we update matrix and Vj according 
to the updated spectral mask matrix W. 



B. P2: capacity optimizing without co-channel interference 

In this section, we assume that the co-channel interference 
is negligible, because the co-channel beams (in frequency do- 
main) or the simultaneously illuminated beams (in time/space 
domain) can be separated far from each other in practice. In 
this way, the capacity optimization can be reduced to a convex 
problem. 

Two cost functions are proposed to solve the frequency and 
time/space capacity optimizing problem without co-channel 
interference. Note that we only discuss the optimization prob- 
lem for BH system because NOFR is dual with BH (see 
Section IV), thus the formulation is also applicable for NOFR 
system by changing the related parameters (e.g. T s — > B c , 
N t — > N c ). 

1 ) n-th Order Difference Cost Function: Here we want to 
match allocated bit rate Ri to requested bit rate Ri as closely 
as possible, i.e., we want to minimize a general function of 
the difference between {Ri} and {Ri} across all the ground 
cells. 



If an n-th order deviation cost function is used, the problem 
can be formulated as 

K 

El ' " 

\Ri — Ri , 

i=l 

subject to Ri,<Ri (20) 

K 
i=l 

iV^ ax is the number of cells illuminated simultaneously, which 
is a satellite payload constraint. N- is the number of time-slot 
allocated to ground cell i. 

We assume that the power allocated to each time-slot 
is constant. And hence Ri can be simplified as (assuming 
Gaussian codes): 



ft = ^S t otlog 2 (l+ 7i ). 



(21) 



Since the co-channel interference is assumed to be neg- 
ligible, the optimization problem shown in (20) is convex. 
Therefore, the lagrangian function is given as 



Let dJ Q^i ^ = 0, we can obtain 



(22) 



RiN t 



B to tlog 2 (l +7i) 



N t 



£totlog 2 (l 



, (23) 



where A is the lagrange multiplier and determined from the 
total available time-slot constraint, which can be obtained by 
solving the equation 



K 



1 Y max J 



(24) 



From (23) and (24) we can obtain 



A = n 



E 



RiN t 



Btotlog 2 (l + 7l ) 



N^Nt 



\ 



E 



N t 



Stot log 2 (l + 7i 



(25) 



If we replace A in (23) with (25), the solution will be 

K 



E 



k=i 



R k N t 



B to tlog 2 (l + ik) 



- NSnNt 



Btot log 2 (l +7 t ) 



(26) 



E 



+ 7i) 
l°g2( 1 + 7fc) 



With the the number of time-slot allocated to each ground 
cell (N-), the throughput allocated to each cell (Ri) can be 
calculated with (21). We should note that the solution in (26) 
is independent of the order n (n > 2) in our case, since we 
suppose BPA, no co-channel interference and the same channel 
attenuation factor (c^) for all the ground cell. 
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2) Fairness Cost Function: Another way to match allocated 
capacity Ri to requested capacity B4 is to maximize the ratio 
between them as 



max 



n(5 

»=i ^ 



subject to Ri < Ri, 

K 

EN 1 < N re N t 



(27) 



i=i 



where is the weighting factor that represents the priority 
of each beam. The problem (27) can be easily converted to a 
convex problem by introducing the logarithm in the objective 
function. Thus, the optimization problem is converted to 



max 



f> iog2 (S- 



(28) 



Thus, the lagrangian function is given as 

J(JV?) = - ^ Wi log 2 ( jpj + A ( ^ Nt - N^Nt] . (29) 



Let = 0, then 



Nt 



uoiRiNt 



(30) 



Aln2Stotlog 2 (l + 7i)' 
K 

With given constraint ^ TV* = N™. dx N t , the lagrange 



i=i 



multiplier can be solved as 

K 

Btot log 2 (l+7«) 

\ = i=l 

iVSaxiV t ln2 • 

The solution will be (replace A in (30) with (31)) 



Nt 



coiRiNt 



N*N t 



f-f log 2 (l + 7fc) 



(31) 



(32) 



Therefore, the throughput allocated to each ground cell (Ri) 
can be calculated with (21). 

VI. Technological Gap 

From Section IV-B we can see that the spectral efficiency 
that each technology can provide makes the real difference. 
Therefore, NOFR and BH systems are not completely dual of 
each other. In this section, we will demonstrate the technolog- 
ical gap between NOFR and BH. Note that we only consider 
the forward (FWD) downlink, because the FWD uplink is not 
a big issue since power at the gateway can be greatly increased 
to compensate the attenuation. 

Equivalent isotropically radiated power (EIRP) is defined as 
(in dB) 



EIRP 



OBO- ^repeater ~~ -^antenna + Gtx , 



(33) 



where Output BackOff (OBO) is the ratio of maximum output 
(saturation) power to actual output power, £ rep eater is the 



repeater loss, L a ntenna is the antenna feed loss, and G tx is the 
satellite Tx. antenna gain. With known EIRP, we can obtain 
FWD downlink C/N (in dBHz) and SNR (in dB) as 

C/N = EIRP - L propagation + (G/T) gt - 101og 1() (fc s ), (34) 

SNR = C/N - 10 log 10 (S c ), (35) 

where ^propagation is the propagation loss, (G/T) gt is the 
ground terminal G/T and kg is the Boltzmann constant. Let 



Psat - L 



repeater 



Gtx ^propagation 



101og 10 (fcs) — 101og 10 (-B c ), and let x\ and x 2 be the OBO 
for NOFR and BH systems, respectively. Therefore, (35) can 
be reformulated as 



SNR f =a-x u or SNR* =a-x 2 - 



(36) 



Let the FWD downlink signal to co-channel interference 
SIR be given as y, Therefore, the FWD downlink SINR can 
be formulated as 



SINR 



down 



SIR' 1 + SNR' 1 = y- 1 + 1Q~^\ (37) 



where x can be x\ or x 2 and SNR can be SNRf or SNR 1 
for NOFR or BH system. 

Let the FWD uplink SINR be z, then the FWD whole link 
SINR is given as 



SINR- Q l = SINR- 1 + SINR^ n 



z 1 + y 1 + 10~. 



(38) 



Let the whole FWD link SINR be 7 = SINR toU the 
spectral efficiency in the case of Shannon limit with Gaussian 
coding can be given as 

i] = log 2 (l +7) ~ log 2 ( 7 ) = -log^z" 1 + 10^), 

(39) 

where we make a high SINR approx given as, log 2 (l + 7) — 
log 2 (7). Therefore, the spectral efficiency for NOFR and BH 
system are 



?/ = -log 2 (z 1 +y 1 + 10 S y t ), 
77* = -log 2 (^- 1 +j / - 1 + 10 f ^). 



(40) 



(41) 



Let the technological gap of spectral efficiency between BH 
and NOFR system Arj be given as 



Arj = rf — rf 



lo S2 — : : ^r- ( 42 ) 



z- 1 +JJ- 1 + 

Let z, x\ and x 2 be constant and x\ > x 2 , Ar] will be a 
monotonically increasing function of y. Therefore, the upper 
bound (maximum) of the technological gap Arj will be 



At] n 



A V \, 



log 2 



;10-(- 



-) 



(43) 



l + zl0-(-To^) 

As we indicated before, the uplink is not relevant. Thus we 
can suppose that the uplink SINR z is constant. The result of 
the technological gap is demonstrated in Fig. 8, it is meaningful 
for us to evaluate the performance of NOFR and BH, and to 
predict the technological gap between NOFR and BH systems. 
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VII. Numerical Results 

The objective of the simulation is: Firstly, to evaluate 
the performance of the proposed capacity optimization ap- 
proaches. Secondly, to compare the proposed system design 
with the conventional design, which is regular frequency reuse 
(fn = 7) and uniform power allocation. Thirdly, to obtain the 
technological gap for a realistic implementation. The payload 
parameters are defined in [24]. 

In order to fairly compare the performance with different 
number of beams in the same coverage (e.g. the European 
countries), we assume that the total traffic requirement is the 
same for all the cases. The linear traffic distribution model 
is defined as Rk — k(3; k = 1, 2, • • • , K, j3 is slope of the 
linear function. The following parameters are assumed in the 
simulations: P sat = 4Watt, B tot = 500MHz, N c = 112, each 
cluster contains 7 beams (as shown in Fig. 2), (3 = 8 x 10 6 bps 
for K = 121, and B c = B tot /N c = 4.4643MHz. 

The parameters of power gain (g), spectral efficiency (rj) 
and traffic matching ratio (p) are studied in the simulation, 
which are defined as the following. 

A. Performance Parameters Definition 

1) Power Gain: We compare the amount of total power 
consumption for joint power and bandwidth optimized allo- 
cation with that for uniform power and bandwidth allocation 
when both achieve the same useful throughput using the same 
total bandwidth. We define the power gain g p as 

KP um 



9p 



K 

fe=i 



(44) 



where P un ; denotes the power per-beam of the uniform allo- 
cation scheme. 

2) Spectral Efficiency: The spectral efficiency is defined 
based on the total allocated traffic and total allocated band- 
width as 



fc=i 

K 

E B * 



(45) 



3) Traffic Matching Ratio: In order to describe the satis- 
faction degree of the allocated traffic with respect to the total 
request traffic, the traffic matching ratio is defined here as 



P = 



K 

E 

fc=i 



R k 



K 



(46) 



E 



Rk 



B. Beam Layout and Antenna Model 

We assume a general beam layout model (shown in Fig. 2). 
A fixed-size space is used to generate different number of 
beams, thus, the beamwidth is decreasing as the number of 
beams increases. It means that the larger the number of beams, 




: z: Cluster 
(~) : Beam Coverage (Cell) 



Fig. 2. Beam layout. 




e ooooooooo 



K=49 
K=121 
K=225 

Toial Throughput Request 



10 15 20 



30 



35 



Fig. 3. Convergence speed. 




K [Number of Beams] 

Fig. 4. Power gain (g p ) vs. number of beams (K). 

the narrower the beamwidth. We assume a tapered aperture 
antenna model with 47.14 dBi maximum antenna gain. Then 
the SINR can be calculated in each iteration of the algorithm 
with a given link budget of a typical Ka-Band (19.95 GHz) 
satellite payload. 

C. Simulation Results 

In order to evaluate the relevance of our iterative algorithm 
(shown in Table II), we perform a study of convergence. It can 
be observed from Fig. 3 that the algorithm is convergent for 
different number of beams, and the convergence is faster with 
the number of beams increasing, e.g. our algorithm runs 24 
and 33 iterations for number of 225 and 49 beams respectively. 



x 




O Shannon case (P1 : capacity optimizing with co-channel interference) 
□ DVB-S2 case (P1 : capacity optimizing with co-channel interference) 

■ O' ' Shannon case (uniform resource allocation) 

■ Q" DVB-S2 case {uniform resource allocation) 



20 40 60 



100 120 140 160 180 200 220 240 260 280 300 
K [Number of Beams] 



1.CQ- 



O Shannon case (P1 : capacity optimizing with co-channel interference) 
— S — DVB-S2 case (P1 : capacity optimizing with co-channel interference) 
O' ' Shannon case (uniform resource allocation) 
'□' DVB-S2 case (uniform resource allocation) 
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Fig. 5. Spectral efficiency (rj) vs. number of beams (K). 



Fig. 7. 



slope (|i)[x 10° bps] 

Traffic matching ration (p) vs. traffic slope (J3). 




1 1 1 1 — 

20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 
K [Number of beams] 



Fig. 6. Traffic matching ration (p) vs. number of beams (K). 

The reason is that the algorithm allocates resources to all 
unsatisfied beams in each iteration, thus, more traffic will be 
allocated with larger number of beams. Consequently, all the 
beams will reach the traffic requirement faster. The algorithm 
has been applied in the realistic satellite payload model and 
proved in [24] that it is applicable to the current multibeam 
satellite system. 

The power gain with respect to the number of beams is 
shown in Fig. 4. We can see that about 6dB and 3.5dB power 
gain can be achieved by capacity optimizing of PI with 
Gaussian coding and DVB-S2 ModCods (defined in [19]), 
respectively (when K — 200). By optimizing the capacity 
achieved per-beam, we do not only reduce the power and band- 
width consumption of small traffic requirement beams, but also 
achieve reasonable proportional fairness from the viewpoint 
of user beams. In Fig. 5, the result shows that the spectral 
efficiency decreases with the number of beams increasing, 
especially when K > 200. The reason is that co-channel 
interference will increase with the beamwidth decreasing. 

In Fig. 6 we can observe that the traffic matches better 
in case of larger number of beams. However, the power 
consumption and the complexity will increase with larger 
number of beams. Therefore, we should balance the total 
achieved throughput with respect to both power consumption 
and complexity. Fig. 7 shows the traffic matching ratio with 
respect to different traffic distribution slope. Obviously, the 
traffic matching ratio drops down with the slope increasing. 
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Fig. 9. Comparison of cost functions in terms of throughput. 

Because the traffic distribution is more asymmetric with larger 
slope. We can see that the capacity optimization can achieve 
better matching ratio compared to the conventional design for 
both Shannon and DVB-S2 cases. 

In order to evaluate the technology gap, we define the 
difference of OBO between NOFR and BH systems as 
Aobo = X\ — x-2- Fig. 8 shows Aobo with respect to A7? max , 
which is defined in (43). We can see that A?7 max is almost 
linear with Aobo, and the slope is increasing with BH system 
OBO (xi) increasing. This result is very useful to predict the 
technological gap between NOFR and BH systems. 

Fig. 9 shows the distribution of throughput for n-order dif- 
ference cost function and fairness cost function along K — 50 
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beams that have a linear distribution traffic demand. In this 
simulation, we assume that (3 = 3 x 10 7 , n = 2 (second 
order function), N t = 32, = 8, the SINR 7fe and the 

weighting factor uik are constant for all the cells in order to 
simplify. The result shows that two different cost functions 
distribute the total available resource (carriers or time/space) 
to all the ground cells with different pattern. Fairness cost 
function is more favorable for low traffic requirement cells 
while n-order cost function distribute more resource to high 
traffic requirement beam. The performance of BH is slightly 
better than NOFR, especially for the low traffic requirement 
beams. Further, the n-order simply neglect too low-loaded 
beams. This is relevant result since it is already considered 
in satcom design. 

VIII. Conclusions 

Current designs of broadband satellite systems are lack of 
the necessary flexibility to match realistic asymmetric traffic 
distributions. Two new technologies are studied to replace the 
current ones over multibeam satellite systems. We prove that 
the two technologies are dual of each other in frequency and 
time domains. Moreover, the technological gap between NOFR 
and BH systems is formulated. Two novel capacity optimiza- 
tion problems, PI and P2, are investigated to best match the 
individual SINR constraints. The current state-of-the art PHY 
layer technology: DVB-S2 and Shannon are implemented in 
order to evaluate the performance. The results show significant 
improvements in terms of power gain, spectral efficiency and 
traffic matching ratio compared to the conventional system. 
For a DVB-S2 and K = 200 case, we can achieve up to 3 dB 
power gain, 0.7 bit/s/Hz spectral efficiency gain, and improve 
10% traffic matching ratio by the proposed capacity optimizing 
approach. For the duality study, the results show that the 
technological gap is only related to the OBO of NOFR and 
BH, and the gap is almost linear with Aobo- Further, we solve 
the second problem P2 with different cost functions. Fairness 
cost function is more favorable for low traffic requirement 
cells while n-order cost function distribute more resource to 
high traffic requirement beams. The study of the capacity 
optimization shows that the BH system performs only slightly 
better than NOFR. We also prove the primary goal of the 
study, that the NOFR and BH technologies can match much 
better than the conventional design in the asymmetrical traffic 
distribution model. 
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